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Abstract — We study the problem of super-resolving a 
superposition of point sources from noisy low-pass data 
with a cut-off frequency f c . Solving a tractable convex 
program is shown to locate the elements of the support 
with high precision as long as they are separated by 2// c 
and the noise level is small with respect to the amplitude 
of the signal. 

I. Introduction 

The problem of super-resolution is of great im- 
portance in applications where the measuring process 
imposes a physical limit on the resolution of the 
available measurements. It is often the case that the 
signal of interest is well modeled as a superposition of 
point sources. Motivated by this, we consider a signal 
of the form 

(U) 



2^ a A 



consisting of a train of Dirac measures with complex 
amplitudes aj located at different locations {tj} in the 
unit interval. Our aim is to estimate x from the lower 
end of its spectrum in the form of n = 2f c + 1 Fourier 
series coefficients ( f c is an integer) perturbed by noise, 



y(k) = 



x(dt) + z(k) 



(1.2) 



for fc e Z, |fe| < f c . To ease notation, we write dl.21 > 
as y = J- n x + z. We model the perturbation z 6 C" 
as having bounded l<i norm, 



< 8. 



(1.3) 



The noise is otherwise arbitrary and can be adversarial. 

Even if we assume the signal x to be very sparse, 
without further conditions to ensure that the support of 
x is not too clustered the super-resolution problem is 
hopelessly ill-posed. This can be checked numerically, 
but also formalized thanks to the the seminal work of 
Slepian J4] on discrete prolate spheroidal sequences 
(see Section 3.2 of []3])- To avoid such extreme ill- 
posedness, we impose a lower bound on the minimum 
separation between the elements of the support of the 
signal. 

Definition 1.1 (Minimum separation): Let T be the 
circle obtained by identifying the endpoints on [0, 1]. 



For a family of points T C T, the minimum separation 
is defined as the closest distance between any two 
elements from T, 



A(T)= inf \t-t'\. (1.4) 



To recover x we propose minimizing the total vari- 
ation of the estimate, a continuous analog to the l\ 
norm for discrete signals (see Appendix A in |3] for 
a rigorous definition), subject to data constraints: 



mm x 



TV 



subject to \\J- n x — y\\ 2 < S, (1.5) 



where the minimization is carried out over the set 
of all finite complex measures x supported on [0,1]. 
For details on how to solve dl.5b using semidefinite 
programming see 0. 

Previous work established that if 



A(T) > — := 2A C 

Jc 



(1.6) 



TV-norm minimization achieves exact recovery in a 
noiseless setting J5]. Additionally, J2) characterized 
the reconstruction error for noisy measurements as 
the target resolution increases. In this work we study 
support detection using this method. If the original 
signal contains a spike of a certain amplitude we could 
ask: How accurately can we recover the position of the 
spike? How does the accuracy depend on the noise 
level, the amplitude of the spike and the amplitude 
of the signal at other locations? Our main result 
quantifies the local stability of the algorithm in order 
to answer such questions. 

Theorem 1.2: Consider the noise model ( 11.31 ) and 
assume the support T satisfies the minimum-separation 
condition ( II. 61 ). The solution to problem ( II. 5I 1 ! 

4 = E &k6 h 

t k et 



1 This solution can be shown to be an atomic measure with discrete 
support under very general conditions. 



with support T obeys the properties 



(i): 



E 



ai 



{t,ef : |£i-tj|<cA c } 

(ii): E 



< C 2 S Vti e T, 



Lemma 2.2: Suppose T obeys condition ( 11.61 a«t/ 
/c > 10. Then for any tj G T there exists a low-pass 



polynomial q tj (t) = Y, k =-f c b k e 



i2 * kt such that 



di | min (tj - ti) < C 3 A C J, 



|i -*,(*)! < 



c[ (t - tjY 



M*)l< 



\ 2 



{tier-. \ii-tj\<c\ c ,tjeT] 
(Hi): ^ \ai\<Ci5, 

{t,ef: |iVt i |>cA c Vt i ST} 

where C\, C 2 and C3 are positive numerical constants 
and c = 0.1649. 

Properties (i) and (ii) show that the estimate clusters 
tightly around each element of the signal, whereas 
(iii) ensures that any spurious spikes detected by the 
algorithm have small amplitude. These bounds are where < c C 2 < C x < 1. 



tj I < cA c , 



(II. 1) 



A? 



\t-ti\ <cA c , U €T/{tj}, 
(11.2) 

\<lt ] {t)\<C 2 , \t-tj\>cX c , tjGT, (113) 



Lemma 12.21 provides a low-frequency polynomial 



essentially optimal for the case of adversarial noise, 

which can be highly concentrated. The main novelty % such that j T q tj (i)x(dt) 
of the result is that it provides local stability guarantees Theorem and the Cauchy-Schwarz inequality 
that only depend on the value of the signal at a 
given location. This allows to quantify the accuracy 
of support detection for each individual spike in the 
signal. 

Corollary 1.3: Under the conditions of Theorem U .21 
for any element ti in the support of x such that ai > 

C26 there exists an element ii in the support of the By Lemmas 12^21 and |2T 
estimate x satisfying 



q t ft)x{dt)~ / q t At)x(dt) 



< Wlu IL \\F n {x ~ x)\\ 2 < 25. 



By Parseval's 



fc 

=-fc 

(H.4) 



U < 



C 3 6 



-C,8 



Despite the aliasing effect of the low-pass filter 
applied to the signal, the bound on the support- 
detection error only depends on the amplitude of the 
corresponding spike (and the noise level). In contrast, 
the local bounds in JE] are not completely decoupled 
from the value of the signal at other locations, which 
makes it challenging to derive similar guarantees. A 
recent paper bounds the support-detection error of a 



Y d k qt j (ik)+ ^ d k (q tj (i k ) 

{fc: \i k -tj\>c\ c } {fc: \i k -tj\<cX c } 

< E \ak\\qt } (tk)\ 

{fc: \t k -tj\>cX c } 

+ E h|i-<mM 

{fc: |ife— tj |<cA c } 

, C[d(i k ,T] ' 



< 



\d k \ min t C 2 , 
i k £T I 



A? 



< CS, (11.5) 



related convex program in the presence of stochastic for a p 0sit i ve numerical constant C. Finally, combining 
noise, but the bound depends on the amplitude of the LemmaEU the triangle inequality, flQ) and we 
solution rather than on the amplitude of the original jj ave 
spike JT]. 

ij - ^ at 

{fc: \i k -tj\<cX c } 



II. Proof of Theorem 1 1.2 1 

We begin with an intermediate result proved in 
Section IH-AI 

Lemma 2.1: Under the assumptions of Theorem U. 2 



^ |afc|min< C a 
t k ef L 



C b d(t k ,T) 



< 25, 



q tj (t)x(dt)- 
E 

{fc: \t k -tj\>cX c } 



q tj (t)x(dt) 
d k q tj (4) 



where C a and Cb are positive numerical constants and 
d (t, T) := min (t — ti) 2 . 

Properties (ii) and (iii) are direct corollaries of 
Lemma ( 12. U . To establish property (i) we need an 
auxiliary lemma, proved in Section III-BI 



+ E a k (q t] {i k ) 

{fc: \i k -t 3 \<cX c } 

< C'5, 

for a positive numerical constant C . 



2 



A. Proof of Lemma 12.71 

The proof relies on a low-pass polynomial provided 
by Proposition 2.1 and Lemma 2.5 in 0. 

Lemma 2.3: Suppose T obeys ( 11.61 l and take any 
v G C' T such that each entry Vj has unit magni- 
tude. Then there exists a low-pass polynomial q(t) = 
Efc=-/ C d k e l2lTkt such that: 

q(tj) = vj, tj e T, 
\q(t)\ <1-C a , 



\t-t j \>c\ c Vtj G T, 



|9(*)I<1- 



A 2 



with < c 2 C b < C a < 1. 

We set vj = ~aj/ \a,j\. The lemma implies 



/ q(t)x(dt)<J2 



a k \ \q(ik)\ 



< J2 1 - min \ 



C b d(t k ,T) 



A? 



\a k 



(II.6) 



The same argument used to prove ( 111.41 yields 



< 25. 



q{t)x{dt) - / q(t)x(dt) 



Now, taking into account that J T q(i)x (dt) = 
by construction and ||i|| TV < 1 1 cc 1 1 T v , we have 



/ q(t)x(dt) 
Jt 



>IN 

> \\x\ 



q{t)x (dt) 

q{t)x{dt) 
,, - 26 



TV 



-2S = J2 



q(t)x{dt) 

28. (II.7) 



«/,■ 



Combining dll.6l ) and dll.7l ) completes the proof. 

B. Proof of Lemma \2.2\ 

The proof of Lemma l2~2l is similar to that of Lemma 



(see Section 2 of 131), where a low-frequency 
kernel and its derivative are used to interpolate an arbi- 
trary sign pattern on a support satisfying the minimum- 
distance condition. More precisely, we set 

q tj {t)= J2 a kK(t-t k )+/3 k KW(t-t k ), (II.8) 

where a, ft G C' T are coefficient vectors, K is 
"sin + l) 7rt)] 4 

K(t)= vv y 1 , teT\{0}, 

(f + l) sin(Trt) 

(H.9) 

and K(0) ~ 1; here, K^ 1 is the £th derivative of 
A'. Note that K, and, consequently, q t are 

trigonometric polynomials of the required degree. 



In order to satisfy the required properties, we impose 

ft,(*i)=l, «t J (*j) = > tj£T/{tj}, 
q' t .(t k )=0, t k e 



T. 



These constraints can be expressed in matrix form. 
Let etj G R' T denote the one-sparse vector with 
one nonzero entry at the position corresponding to tj . 
Then, 



~D 


Di 




a 




V 




D 2 




A 








r tj \ < cA c , ^ G T, where 



(D ) lk = K (ti — t k ) , {D{) lk = (< ; - t fc ) , 
(^2) ;fe -A( 2 ) (t,-i fc ), 

and 2 and fc range from 1 to \T\. It is shown in 
Section 2.3.1 of [31 that under the minimum-separation 
condition this system is invertible, so that a and j3 
are well defined and q(tj) = 1. These vectors can be 
expressed as 



S- 1 ^, S:=D -D 1 D^ 1 D 1 , 



where S is the Schur complement. Let HMHoo denote 
the usual infinity norm of a matrix M defined as 

\\M\loo = maxi^H^! HMxIloo = max; £] fc \M lk \. 
Borrowing some results from Section 2.3.1 in ||3], 



a 




1 


A 




r D 2 l Di_ 



I'-Slloo < 8.74710" 3 , 
I IS 1 " 1 1 1 < 1 + 8. 824 10" 



lll-S*" 1 !! <||S , ~~ 1 || \\S - 1 

II I I OO — II IIOQ 11 

lla-edl <||/-S'" 1 || 

II l 3 I I OO — I I II 

< 8.825 10~ 3 , 
Halloo < 3.294 10~ 2 A C . 

Lemma 2.6 in allows to obtain 



< 8.82510" 



e tj\ 



(11.10) 
(11.11) 



1 47T 

K (t) < — — j < 0.333, K (t)< — < 4.18 f c 
for \t\ > cA c as long f c > 10. By the same lemma, 

E l# 

1 



OO 

^E 



;= - (/ c A mm (l/2 + 0) 4 

E 



E 



1 



(/c^minO 



< 1.083, 



tk£T\{t a ,t b } 



4?r 



<y — 

-t^/ c 3 (A mm (l/2 + /)) 4 



4tt 



E 

j_0 /c (^minO 



< 1-75 f c 



where t a and % are the two spikes nearest to t. Let 
ti be the element of Tj {tj} that is nearest to t. 



3 



Combining these inequalities with ( III. 101 ) and ( III. 1 II ) 
proves that 

= | E akK (* - + E /3fei ^ (1) (< - 

t fc 6T t fc GT 

t fc eT\{ti ,*.,■} 

+ Plloo(|^ (1) (i-t^l + l^ 1 ' (t-ti) 

+ (*-**) 

t k eT\{ti,tj} 
< 0.69, if \t-t k \ > c\ c , t k g T. 



so that (III. 3b holds. The proof is completed by two 
lemmas which prove (III. lb and ( III. 2b and \q tj (t)| < 1 
for any t. They rely on the following bounds borrowed 
from Section 2 of (3), 



K{t) > 0.9539, 



K {2) (t) < -2.923 / c 2 



if« (t) 
if( 3 ) (t) 



< 0.5595 f c , 

< 5.697 / 3 , 



if < 2 > (t) 



< 3.393 / 2 , 



(11.12) 

and on the fact that, due to Lemma 2.7 in J3) it holds 
that for any to g T and t e T obeying |t — io| < cA c , 



E |^ (2) (*-m 

tfeeT\{t } 

E \K {3 Ht-h) 

t k eT\{t„} 



<1.06/ c 2 
< 18.6/ 3 . 



(11.13) 
(11.14) 



Lemma 2.4: For any t such that \t — tj \ < cA c , 

, 4.07 (i-i,-) 2 , 2.30 (t-f,) 2 
1 ^ ^- < ft (f ) < 1 i — 



A, 2 



A, 2 



Proof We assume without loss of generality that tj = 
0. By symmetry, it suffices to show the claim for t g 
(0, c A c ]. By dTTTOb . dlLTTi dTTTZb . flLOb and (ITUl) . 

qg (t) = E "fe^ (2) (* - **) + E ^ Ri3) (* - **) 

tfcST t fc £T 

< (l - I la - e t . 



+ ||"- % 

+ ll/3|lc 



I loo E l^-M 

t fc eT\{o} 

i^(t)| + e |# (3) 

i t er\{o} 
< -2.30 / 2 . 

Similar computations yield \q$ (t)| < 4.07 J 2 . Since 
9o (0) = 1 and 9o (0) = 0, this implies 1-4.07 {f c t) 2 < 
q Q (t) < 1-2.30 (f c t) 2 . * 



Lemma 2.5: For any ti g Tj {tj} and t obeying 
\t — tjl < cA c , we have 



M*)l< 



16.64 (t-tj)' 
A? 



Proof We assume without loss of generality that t; = 
0. By symmetry, it suffices to show the claim for t g 
(0,cA c ]. By flrTOl flLTTT ). ( llTT2b . flU3b and flLTl . 
we obtain 



E^ (2) (t- tk )+j2^K {3) (t-t k ) 



t k £T 



t k £T 



< O + lh+^IIJ ^ (2) 



a — e+ . 



x( 2 ) (t) 



E |* (2) (*-**) 

t k er\{o} 



(t) 



E |^ (3) (*-4) 

t k er\{o} 



< 16.64 / c 2 



where we have also used that in the interval of interest 
\ K(2) (t-t 3 )\ < ^(A^ie/^ < 15-67 / 2 due to 
Lemma 2.6 in 0. Since qo(0) = and q' Q (0) = are 
equal to zero, this implies |(?o(t)| < 16.64(/ c t) 2 . ■ 

Acknowledgements 

This work was supported by a Fundacion Caja 
Madrid Fellowship. The author is grateful to Em- 
manuel Candes for useful comments regarding this 
manuscript and for his support. 

References 

[1] J. M. Azais, Y. de Castro and F. Gamboa. Spike detection from 

innacurate samplings. Preprint. 
[2] E. J. Candes and C. Femandez-Granda. Super-resolution from 

noisy data. Preprint. 
[3] E. J. Candes and C. Fernandez-Granda. Towards a mathematical 

theory of super-resolution. Communications on Pure and 

Applied Mathematics, to appear. 
[4] D. Slepian. Prolate spheroidal wave functions, Fourier analysis, 

and uncertainty. V - The discrete case. Bell System Technical 

Journal, 57:1371-1430, 1978. 



4 



